Semiclassical Quantization for the Spherically Symmetric Systems under an 

Aharonov-Bohm magnètic flux 
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The semiclassical quantization rule is derived for a system with a spherically symmetric potential 
V(r) ~ r" (—2 < v < oo) and an Aharonov-Bohm magnètic flux. Numerical results are presented 
and compared with known results for models with u = — 1,0,2, oo. Itis shown that the results 
provided by our method are in good agreement with previous results. One expects that the semi- 
classical quantization rule shown in this paper will provide a good approximation for all principle 
quantum number even the rule is derived in the large principal quantum number limit n ^> 1. We 
also discuss the power parameter v dependence of the energy spectra pattern in this paper. 
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I. INTRODUCTION 



In the past 20 years, Aharonov-Bohm (A-B) effect, a topological non-local physical effect at the quantum level, has 
been of much interest in the studies of còsmic string [1], (2 + 1)-D gravity theories [2] and especially in the context 
■^j- of anyon [3], which has shed light on the understanding of the phenomenon of the fractional quantum Hall effect 
. [4-7], superconductivity [7,8], repulsive Bose gases [9], and so forth. There are only a few models coupled to different 
potentials along with an A-B magnètic flux that can be solved exactly. For the system with both an A-B magnètic flux 
and a spherically symmetric potential of the form V(r) = Ar' y (— 2 < v < oo), the solvable models known to us include 
the cases with the parameter v = —1,0, 2, oo [10-12, 16]. Here A is a constant parameter. Note that when v = — 1, it 
is a system with both an A-B magnètic flux and a Coulomb potential (A-B-C) [10-12]. This system describes the 
■-«^ • relative motion of two charged partides, with one of them carrying elèctric charge and magnètic flux (—q,—<ï>/Z) 
i-S^ | while the other one carrying (Zq, <£>). Here Z{^ 0) is a non-vanishing real number. This system is of much interest in 
ípH. many different areas [12]. 

In the past three dècades, much progress has been made in the semiclassical methods toward the understanding of 
these systems. These kinds of semiclassical methods provide with us a powerful approximation tool in different areas in 
order to extract useful information from various unsolved problems including the quantization of the classical chaotic 
systems [13], deformed atòmic nuclei, asymmetric fission nuclei [14], semiclassical quantum dots, and weak localization 
in mesoscopic systems [15]. In this paper, we will consider a generalized system with both an A-B magnètic flux and 
a spherically symmetric potential of the form mentioned above. The set of the parameters (A, v) will be discussed in 
the following ranges (i) (A < 0, —2 < v < 0) and (ii) (A > 0, v > 0). We will derive a semiclassical quantization rule 
of the approximated energy spectra for this set of parameters. The distribution tendeney of the energy spectra on 
different vàlues of the parameter v will also be given. By comparing with the known results, including the models with 
v = —1,0,2, we find that our method agrees with these exact results. In addition, for the exactly solvable model with 
v = oo, the difference between the exact and semiclassical results will be shown to be very small from a numerical 
computation. Therefore, we are confident in that our formulae will also provide a good approximation for the two 
ranges of parameters mentioned above where v ^ — 1, 0, 2, oo. 

This paper is organized as follows. In section II, we will derive the semiclassical quantization rule of the A-B 
effect under a spherically symmetric potential. In particular, we will first derive the non-integrable phase factor of 
the Green's function due to the A-B effect in a spherically symmetric system. The corresponding radial Schròdinger 
equation will also be derived accordingly. The semiclassical wave functions will also be derived according to the 
semiclassical consideration of the Bohr's corresponding principle. Consequently, the quantization rule can thus be 
obtained by comparing with the well-known WKB phase. We will also study the distribution dependence of energy 
spectra in various models in section III. The effect of magnètic flux will also be discussed and emphasized in this 
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section. Finally, in section IV, some conclusions will be drawn. In order to provide a self-contained information, wc 
will show the WKB matching condition of the semiclassical wave functions in the appendix. 



II. SEMICLASSICAL QUANTIZATION RULE OF THE A-B EFFECT WITH A SPHERICALLY 

SYMMETRIC POTENTIAL 



The fixed-energy Green's function G°(r, r'; E) for a charged particle with mass m propagating from r to r' satisfics 
the Schròdingcr cquation 



E - H (r, -V) 



G°(r,r';E) = <5 3 (r-r'), 



(1) 



where the system Hamiltonian is given by H = — h 2 V 2 /2m + V(r) as usual. In the spherically symmetric cases, the 
angular dccomposition of the Green's function can be written as 



oo / 



G°(r,r'-E) = J2 £ E)Y lk (6, ,<f') 



(2) 



1=0 k=- 



with Yik the well-known spherical harmonies. As a result, the lcft hand side of Eq. (1) can be brought to the following 
form 



oo l 



1=0 k=-l 



d 2 2 d \ t 1(1 + í)h 2 



2m \ dr 2 r dr 



2mr 2 



V(r) 



xG»(r/;£)y lt (^K(9>'). (3) 
For a charged particle in a magnètic ficld, the Green's function G is related to G° by the following equation 

G(r,r';E) = G°(r,r';E)e£iï A{í) - dí , (4) 

with a globally path-dependent non-integrable phase factor [17,18] given above. Here we have used the vector potential 
A(r) to represent the magnètic field. For the Aharonov-Bohm magnètic flux under consideration, the vector potential 
can be written as 

where the two-dimensional radial lcngth is defined as p 2 — x 2 +y as usual. Moreover, è v is the unit vector of coordinate 
ip, e is the radius of region where magnètic field exists. Hence the total magnètic flux is given by $ = ire 2 B. Note 
that the associated magnètic field lines are confined inside a tube, with radius e, along the z-axis. Along the region 
without magnètic field, the path-dependent non-integrable phase factor is given by 



-ipo d\' tp(X') 



(6) 



where we have used the subscript P to represent the path dependent nature of phase factor and we have denoted 
(p(X') = dip/dX' . Also, /j,q = —2eg/hc is a dimensionless number defined by <ï> = 47rg. The minus sign is a matter of 
convention. According to the discussion in Ref. [18], only phase factors with closed-loop contour are considered where 
the description of electromagnètic phenomenon are complete. Hence, we have 



i- £d\'iï\>) (7) 



with integer vàlues n corresponding the winding number. The magnètic interaction is therefore purely topological. 
Therefore the nonintegrable phase factor becomes 

e -iHo(2nir) ^ ^ 
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With the help of the equality between the associated Legendre polynomial PJf(z) and the Jacobi function P^^\z) 
[19,20], we find that 



Therefore the angular part of the Green's function in the expression (3) can be turned into the following form 



(9) 



E 

fe=-i 



2/ + 1 r (z - + 1) r (z + fc + 1) 



47T 



r 2 (z + i) 



e e 1 e e 1 

cos - cos — sm - sm — 
2 2 2 2 



x p/_\ fe) (cos 9)P^ k k) (cos fe (^') . 



(10) 



In ordcr to include the non- integrable phase factor due to the A-B effect, we will change the index l into q related by 
the definition l — k = q. As a result one can rewrite the Eq. (3) as 



oc oc 

*-E E 

K q=0 k= — oo 



2m V dr 2 r dr 



{q + k){q + k + l)h 2 
2mr 2 



V(r) 



xG° q+k (r/;E) 



2(q + k) + ÍT (q + Ï)T (q + 2k + ï) 



47T 



r 2 (q + k + 1) 



cos - cos — sm - sm — 
2 2 2 2 



x P^ (cos 6)P^ (cos 0')e ífc (^') . 



(11) 



In addition, the non-integrable phase in Eq. (8) can now be included with the help of the Poisson's summation formula 
(p.124, [21]) 



oc «oc oo 

f(k) = dy ^ nyi f{y). 

_ „„ J — OO „ _ w 



/e— — oo 



(12) 



Therefore, the expression (11) can be written as 



q—O k— — oo 



2m V dr 2 ^ r dr 



2 d \ (q + z)(q + z + í)h 2 



2mr 2 



-V(r) 



xG q+z (r,r';E) 



2(q + z) + ÍT (q + 1)T (q + 2z + 1) 



4tt 



T 2 (q + z + í) 



9 9' 9 9' 
cos - cos — sm - sm — 
2 2 2 2 



x P^ (cos 9)P^ (cos 0>y(>-K>)(v+2k*- v ') ^ 



(13) 



where the superscript in G q+k has been suppressed to reflect the inclusion of the A-B effect. The summation over 
all indices k forces z = /Uq modulo an arbitrary integer numbcr. Therefore, one has 



oo oo 



*-E E 

g=0fc=-oo 



(<P_ 2_d\ (g+|fc + Mo |)( g +|fc + /xo| + l)fi 2 
2m \ dr 2 r dr j 2mr 2 



V(r) 
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, r , ,. p s í [2 (g + |fc + Mol) + 1] r (g + 1) r (2 |fc + Mo | + g + 1) ) Jk ( v _ v >) 
.G g+lk+ , ol (r,r,E)^ p2 - - - - - - je l J 



x (cos0/2cos0y2sin0/2sin072) |fe+ ' ío1 P^+^l^+^D^os^P^+^l·l^l^cos^). (14) 

Note that the influence of the A-B cffcct to the radial Green's function is to replace the integer quantum number l 
with a fractional quantum number q + \k + Mol- Analogously the same procedure can be applied to the delta function 
S 3 (r — r') in the rhs of the Eq. (1) with the help of the following sòlid angle representation of the 6 function 



oo l 



í(íi-n') = E E ifeOMWV)- 



(15) 



Z=0 k=- 



Therefore, for the set of the fixed quantum numbers (q, k) one can show that the radial Green's function satisfics 



E 



— + -—\ + (g + l fc + HXg + |fc + Mol + i)fi 2 

2m V c^ 2 r cír / 2mr 2 



V(r) 



x G q+lk+flol (r,r';E) = 5(r-r'). 
As a result, the corresponding radial wave equation reads 

H 2 7 ( 7 + l) 



(16) 



2m dr 2 



E - 



V(r) + 



u 7 (r) = 0, 



2m r ; 

where we have set 7 = q + \k + Mo| 5 and w 7 (r) = ri?^ 7 (r). Obviously, i?s 7 satisfies the spherical Bessel equation 

7(7+1) 



— + -— + I ( 1 

dr 2 r dr \ r 2 



Rnj(r) = 



(17) 



(18) 



with the definitions k = y2mE/h 2 and the reduced potential U(r) = 2mV(r)/h 2 . For simplicity, we have written 

Rn 1 {r) instead of Rü, q ,k(r) in which each set (n,q,k) denotes a quantum state. Hence the A-B effect reflects itself 
by the coupling to the angular momentum in radial Green's function which turns the integer quantum number into 
a fractional one. 

To find the semiclassical quantization rulc, lct us first consider the asymptotic form of the bound state wave 
functions of a charge article moving in a spherically symmetric potential of the form V(r) = Xr l/ (X < 0, — 2 < v < 0) 
undcr an A-B magnètic flux for the energy limit E — > 0. Due to the Bohr corresponding principle, this stands for the 
semiclassical approximation since there are infinitely dcnscd cncrgy levels near E — > 0~. According to the Eq. (17), 
the asymptotic wave equation reads, in the E — > CP limit, 



We can also perform the following transformations: 



ft 2 7(7+1) 
2m r 2 



w 7 (r) = 0. 



\ ,u(r)=W(p). 



Consequently, the Eq. (19) yields 



d 2 W 
Ip 2 ' 



7(7 + 1) 



which can be further reduced with the help of the following changc of variables 

2 



i/ + 2' 



>+2)/2 



W(p) 



,1/(^+2) 



v(z). 



(19) 



(20) 



(21) 



(22) 
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As a result, the Eq. (21) becomes 



d v 1 dv 

dz 2 z dz 



2 7 +l 
u + 2 



v = 0. 



(23) 



This is exactly the Bessel's equation of integral order vi = (2-f + 1) / {y + 2). The boundary condition (B.C.) of the 
function u(r) in the Eq. (17) is simply u(0) = 0. Therefore, the corresponding B.C. of v(z) in the Eq. (23) is v(0) = 0. 
The Bessel function of the first kind is known to be the solution of the Bessel's equation. Therefore, by imposing the 
B.C. appropriately, one can show that 



v(z) = J Vl (z) 



(24) 



is the solution of the Eq. (23) with the prescribed boundary condition. Therefore, the solution of the radial wave 
equation near E — > becomes 



í (r)=W(p)=z 1 /^J Vl (z). 



(25) 



^From the asymptotic behavior of the Bessel function near r — > 0, or equivalently p — > and z — > 0, one can show 
that 



u(r) 



z l/{u+2) z (2 1+ l)/(u+2) _ z (2~t+2)/(v+2) _ ^+1 _ 



On the other hand, from the asymptotic behavior of the Bessel function approaching r — > oo, 

T 

one can show that 



/ 2 / af f\ 



(26) 



(27) 



-^ /4 -(^ +2)/2 --f- 4 )- ^ 

Note that one can also compute the following integral, near the limit E — > 0, and show that the following identitics 
hold: 

1/2 



h 2 J v + 2' z/ + 2 

It follows that, in the limit E — > and r — > oo 



\ 1 /2 

2m l A l\ 2 r (^+2)/2 _ 2 p (v+2)/2_ 



u(r) <~ r cos 



~ r sin 




2m ,_ , .. 7T 7T 



(30) 



where V(r) = Xr u (X < 0,— 2 < v < 0). If we take the integration upper bound r as the classical turning point r c 
where V(r c ) — E , the phasc of u(r) can be shown to be the WKB phase (see Appendix A for details) 



u(r c ) oc sin 



n + 



(31) 
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Consequently, from comparing the equations (30) and (31), one can extract the following quantization condition 



v / 2m{E -V(r)) 



n + 



2 7 + ^ + 3 
2(i/ + 2) 



irK, for n = 0, 1, 2, 3 • 



(32) 



Here n is the radial quantum number. Although the Eq. (30) is obtaincd in the limit E — > 0, or cquivalently in the 
large quantum number where n » 1, above result can still be extended to all possible vàlues of n. In fact, the integral 
in the Eq. (32) can be written in terms of analytic form. Indeed, with the help of the following change of variables 



r v = esc 2 £, 



one can re-write the following integral as 

£^2m{E-V(r)) = ' V^W\ cos 2 Ç(smÇ)-V·'- 2 dÇ. 

In addition, with the help of the following formula (see, for example, Ref. [23], p8), 



one has 



/ 

Jo 



rr c o f E\ 1 t V 

j y/2m(E-V(r)) = -- ^-J V / 2^\E 



2I> + < ? )' 

^ r(-i-j) 
4 r(i-i) ' 



Inserting the result of the Eq. (36) into the Eq. (32), one has 

E -- 



i^ 2/ " +2) (y 



(33) 



(34) 



(35) 



(36) 



2 \u\ ^/tt ( n + 



2(<7 + |fc + /i |) + z/ + 3\ r(l 



2í/ + 4 



-i 2^/(^+2) 



r(-i-è: 



(37) 



where the ranges of the parameters are A, S < 0, — 2 < v < 0, n, q = 0, 1, 2, 
with the potential of the form V(r) ~ — e 2 /r, one has 



E n , q ,k = -mc 2 



a 



2[n + q+ |fc + ^ | + 1] 



2 • 



and — oo < k < oc. For example, 



(38) 



Here a = e 2 /Tic denotes the fine structure constant. This agrees with the exact result given in the Ref. [11]. We 
see that the A-B effect has changed the splitting of energy levels although the electron moves in the absence of the 
magnètic field. In addition, when the flux is quantized, namely, Ang — {2·nhc/e)x integer, \k + fj, \ is an integer and 
hence the spectrum is the same as the energy spectrum of the pure hydrogen atom. 

To obtain the semiclassical quantization rule for all positive power v > of the potential V(r) = \r v , one can 
perform the following change of variable 



p = r a , u(r) — p l3 v(p) 



(39) 



and show that the Eq. (17), becomes 



dr 2 



2 n 2+p-2/ a d 2 v(p) 



a p 



. + a 2 (2p + 1 - I)pi+0-V«^) + a 2 P(P - I)/- 2 /%(p). 
dp z a dp 



a 



(40) 
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Note that the different ranges v > and v < can be proper ly adjusted when the parameters a and (3 are chosen 
appropriately [22]. In addition, if we set 



v 1 (■> v ' 

« = ~7. /3 = -o 1 + - 
f' 2 V f 



(41) 



the term dv/dp in the Eq. (40) disappears. Inserting this back into the Eq. (39) and then the Eq. (17), one has 



2m 



dp 2 



E I - 

v 



2m f 



7(7 + 1) 



4 W 



(42) 



If we choose 1/ = — 2i// (2 + 1/), the above equation reduces to 



2m dp 2 



E' - \'p v ' - i{i + 1) 



2mp 2 



v = 



with the following relations linking different parameters 



( v' = -- 2 - 



(2+f) 

# = "A (£) 



l ' W ~ 2 / l/ 2 i/+2 2 



(43) 



(44) 



Note that the structure of the Eq.s (17) and (43) are similar except the signs of the parameters. Accordingly, the 
eigen solutions for A, v, E > can be found from X',i/',E' < 0. Inserting the relations (44) into the Eq. (37), one 
thus finds that 



E = A 2 ^+ 2 ) f —) 
\2mJ 



\ v /=( n , (g+lfc + Mol) , 3\ r(i + |) 



2v/{v+2) 



(45) 



with the ranges of the parameters A, v, E > 0, n, q = 0, 1, 2, • • • , and —00 < k < 00. As a realization, the three- 
dimcnsional simple harmònic oscillator moving in the presence of the A-B magnètic flux can be described by the 
model with the parameters v = 2 and A = rnuo 2 /2. Hence we can calculate the energy engenvalue from the Eq. (45) 
that gives us the following result 



2n+ (q+ \k + no\) ■ 



Another example is given by the model with an infinitely deep potential 

V(r) = \r» = | ~' 



heu. 



for r > a 
for r < a 



Similarly, Eq. (45) implies the following energy spectra 

h 2 TT 2 f q+ \k + fi \ 



E, 



n,q,k 



2ma 2 



(46) 



(47) 



(48) 
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Here we have replaced (n+7/2+3/4) with (n+7/2+l) according to the matching condition of the WKB approximation 
given in the Appendix A. The analytic energy spectra of this system is then given by the zeros of the modificd Bessel 
function in the Eq. (3.28) of the Ref. [16] 

fy/-2mE \ , s 

I q +\k+no I+1/2 I ^ a 1=0. (49) 
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OC 

UJ 
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FIG. 1. Comparison of the exact and approximate energy eigenvalue as a function of the radial quantum number n. (a)The 
exact and approximate energy eigenvalue are shown in Fig. l(a). (b)Their difference is shown in Fig. l(b). Here we have set 
q + \k + fi \ = 2.5. See Eq.s 48 and 49 for details.The unit of the energy eigenvalue is set as h 2 n 2 /2ma 2 . 



The numerical analysis shown in Fig. 1 (a) (for q+ \k + no\ = 2.5) indicatcs that the result (48) is in good agreement 
with the exact result (49). In addition, Fig. 1 (b) exhibits the difference between the exact and approximate results. 



III. THE 1/-DEPENDENCE OF THE DISTRIBUTION OF THE ENERGY SPECTRA 

Note that the Eq. (45) indicates that 

/ , g +|fc + Mo| , 3^ +2) 
E n ,g,k oc I n+ h - 

For example, for the model with an inhnitely deep potential (i.e. v — > 00), one has 

/ q+\k + hq\ 
E n , q , k oc n + ' + 1 



(50) 



(51) 



On the other hand, from the Eq. (37), one has (when — 2 < v < 0) 
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■^n.q.k 0^ 



n + 



2(q + \k + no\) + + 3 



2u/(v+2) 



2v + 4 

In addition, we can calculate their derivative with respect to n and find that 



dn 



> 



(52) 



(53) 



for all considered models. Thus one expeets that the energy levels E n ^ will monotonically increase as n monotonically 
increases. The q and |fc + ^ | dependence of the energy eigenvalue £7 n ,g,fc can be found by the Hellmann-Feynman 
formula (e.g. [24]) 



dE nq u 



dq 



dH 



dq 



where the Hamiltonian is given by 

t 

2m dr 2 

Thus, we can derive the following results: 

dE rhq ^ k 



U : : + ( xr" + ^ (g+l fc + Mol)(9+|fc + Mol + l) 

2m r 2 



dq 
d \k + Mol 



n,q,k 



[2(g+|A! + /io|) + l] h 2 



q . k 



2mr 2 

\k + no\) + l] h 2 



2mr 2 



n,q,k 



> o. 



(54) 



(55) 



(56) 



(57) 



This means that the energy spectra E n , q ,k will monotonically increase as any one of the quantum numbers in the set 
(n, q, k) monotonically increases. Therefore the ground state will be given by n = q = k = 0. The details can be 
obtained by analyzing the tendeney of E n ^\. with respect to the change of the parameter v. 





FIG. 2. Energy as a function of q and n for four difïerent u's are shown in Fig. 2. Here we have chosen \k + /Uo| = 0.5. The 
unit of the energy eigenvalue is set as mc 2 a 2 /2, (9ir 2 \ 2 h 2 /8rn) 1 3 , fiu), and h 2 iv 2 /2ma 2 for Fig.2(a), Fig.2(b), Fig.2(c), and 
Fig.2(d) respectively. 
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A. Distribution Tendency of The Energy Spectra for v = — 1 



The energy spectra for a charged particle moving in the Coulomb potential and an A-B flux is given by the Eq. 
(38). Its first and second order derivatives with respect to the parameters (n, q, \k + /Uq|) are 



dn 
d 2 E n 
dn' 2 



(58) 



= mc 2 a 2 



a« a - mc a («+ g +|fe+ M0 |+i) 4 < u 



(59) 



ggg^ _ TOC 2 Q! 2 i r > 



8|fe+í*or 



9 9 

mc cr 



(n+rj+lfc+pol+l) 



T < 



(60) 



Consequently, E n ^ q ^ tends to increase and saturate gradually as anyone of the parameters in the set (n,q, \k + /xo|) 
increases. It implies the bending curve as shown in Fig.2 (a). The unit of the energy eigenvalue in Fig. 2 (a) is chosen 
as mc 2 a 2 /2. 



B. Distribution Tendency of the Energy Spectra for v 

The energy levels for the model with v = 1 are given by the Eq. (45) 



A 2 7^ 1/3 



2m 



3tt / ( g +|fc + M |) , 3 
~2~ l 2 + 4 



2/3 



Their derivatives with respect to the parameters (n,q,\k + fio\) yield 



(61) 



On 
d^ 2 



7T 

1/3 



3jr / , (g+|fc+Mol) , 3 
2 ' T 2 ^4 



-1/3 



>o, 



^ 2m ) \2 ) [2 \ 



-4/3 



< 



(62) 



dE„ 



dq 



2m 



1/3 



2 

1/3 



(q+|fc+M |) 
2 



¥( 

(«)"*(*) [*(" + 



-1/3 



> o, 



(q+|fc+£o|) 
2 



-4/3 



< 



(63) 



and 



dE niqtk 
d\k+fi a \ 

d 2 E n ^ k 



/ 2 2 \ I/ 3 

' ■? 9 \ 1/3 
2m V 8 



( g +|fc+ M o|) 



+ 



-1/3 



3_7T 

2 



(q+|fc+Mo|) 
2 



>o, 

-4/3 



< 



(64) 



It is obvious that E n ^ q ^ will monotonically increase when the value of any parameter of the set (n, q, \k + fio\) increases 
as shown in Fig.2 (b). Note that the slope is much more smooth than the model with v = — 1. The unit of energy in 

Fig.2(b) is chosen as (9ir 2 \ 2 h 2 /8m) 1/3 . 



C. Distribution Tendency of the Energy Spectra for v — 2 



The energy spectra for a charged particle moving in the three-dimensional harmònic potential and an A-B flux is 
given by the Eq. (46). Its first and second order derivatives with respect to the set of parameters (n, q, \k + Ho\) read 
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2hco (const.) 



dn 



dn 2 



(65) 



dE 



n,q,k 



and 



dq 

dE ntqí k 
d\k + /i | 



ha; (const.), 



hui (const.), 



J n,q,k 



8q 2 



= 0, 



d 2 E. 



n,q,k 



d\k + (j,o\ 



(66) 



(67) 



This means that E nq j. will linearly increase as any one of the parameter in the set (n,q,\k + u \) increases. The 
details is shown in Fig.2 (c) with the unit of energy given by fruj. 



D. Distribution Tendency of the Energy Spectra for v = oo 



According to the Eq. (48), we obtain the first and second order derivatives with respect to E n ^k 



dn 



2j;2 r 



ir 2 h 



ma* 



(q+ |fc + /i | 



+ 1 



> 0, 



d 2 E n: g :k 
dn 2 



2*2 



ir 2 h 



ma* 



>0, 



(68) 



dE, 



n,q,k 



Dq 



and 



dR 



n,q,k 



d\k 



Mol 



7r 2 h 2 
2ma 2 



ir 2 h 2 
2ma 2 



n+ (í±J*±M + 1 



n + 



Mol 



+ 1 



>0, 



>o, 



d 2 E 



n,q,k 



dq 2 



d 2 E niqi k 



TT 2 h 2 

Ama 2 



>0, 



2*2 



TT 2 H 



d\k + u \ 2 ^ma 2 



>0, 



(69) 



(70) 



for the model with v — > oo. Note that E n ^k will monotonically increase when any one of the parameters in the set 
(n, q, \ k + no\) increases. The rate of increase is, however, faster than the model v = 2 since the curve climbs up as 
shown in Fig.2 (d) with the unit chosen as h 2 ir 2 /2ma 2 . 

In summary, all thcse results imply the following rules for a charged particle moving in the spherically symmetric 
potential V(r) = \r v (—2 < v < oo) [22] and an A-B magnètic flux: 

(a) The energy spectra of the bound states depend on the quantum numbcr (n, q, fe) and monotonically increase as 
any one of the quantum numbers increases. 

(b) whcn v = 2, the energy spectra E n ^ q ^ depend linearly on any parameter in the set (n,q, fe); when v > 2, the 
energy curve bends up as any one of the quantum numbers (n, q, fe) increases. On the other hand, when v < 2, the 
curve bends down as any one of the quantum numbers increases. 

(c) when v = 2, we h&vedE/dn : dE/dq = 2:1, dE/dn : dE/d\k + u \ =2:1, and dE/dq : dE/d\k + u \ = 1:1 
which are related to the closeness of the classical orbits and whether the model is exactly solvable or not. For the 
case with positive powcr of v , 



E 



n + 



k + Mol) 3 
2 4 



2i//(i/+2) 



Although we still have the same ratio of derivatives, the above relation does not hold for the exact solution. 

(d) whcn v = —1, its energy spectra have the properties, dE/dn : dE/dq = 1:1, dE/dn : dE/d\k + uq\ =1:1, 
and dE/dq : dE/d\k + u \ = 1 : 1. They are also related to the closeness of the classical orbits. For the models with 
negative power of v, (—2 < v < 0), the WKB approximation given by the Eq. (37) implies that 



E 



2(q+ \k + no\) + + 3 
2v + A 



-i 2v/(v+2) 



This hence implies that dE/dn : dE/dq = v + 2, dE/dn : dE/d \k + fj, \ = v + 2, and dE/dq : dE/d \k 
are all equal. This relation does not hold for the exact result for the same reason. 



Mol 
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(e) The increase of intensity of thc magnètic flux will change the slope of thc encrgy distribution in both thc models 
with — 2 < v < and the models with < v < oo. More explicitly, when v < 2, increasing the flux will depress the 
slope; whereas when v > 2, increasing the flux will lead to the increase of the slope. In addition, the model with 
v = 2 is marginal in the sense that the slope of the energy distribution will not be affected by the change of the flux. 
For details, see the difference shown in the Fig.2 and Fig.3. Note that \k + /M)\ is set as 0.5 and 12 in Fig.2 and Fig.3 
respectively. 



(a) v = -1 (b) v = 1 




FIG. 3. Energy as a mnetion of q and n for four different i/s. Here we choose |fc + /io| = 12. The unit of the energy eigenvalue 
is set as mc 2 a 2 /2, (9n 2 \ 2 h 2 /8rn) 1 3 , hu>, and Tl 2 -ïï 2 /2ma 2 for Fig.3(a), Fig.3(b), Fig.3(c), and Fig.3(d) respectively. 



IV. CONCLUSION 



The semiclassical quantization rule is presented for a charged particle moving in a system with a general central 
force described by the potential V(r) = \r v , with — 2 < v < oo, and an A-B magnètic flux. The formulae obtained 
in this paper are in good agreement with the energy levels with all known exactly solvable models with some specific 
vàlues of v. Furthermore, we have presented numerical results for v = oo which is also in good agreement with the 
exact result. Therefore, one expeets that the semiclassical quantization rules will also be in good agreement with the 
models prescribed by a large ranges of v even the results shown in this paper are more reliable for the case with large 
principle quantum number n. 



V. APPENDIX A 



The WKB wave function for a charged particle moving in a smooth potential well near the neighborhood x ~ a 
(x > a), where x = a,b are the intersection points of the horizonal line y — E and the curve y = V(x) as shown in 
Fig. 4(a), can be expressed in terms of the classical momentum p as (see, for example, Ref. [24] for details) 



C 

^>(x) = sin 



pdx + — 



— sma(x), 



(71) 



12 



were C is constant. Analogously, near the neighborhood x ~ b (x < b) we have 



c 

ifi(x) = — — sin 

Vp 



1 í b J 


ir 


— / pdx 4 


' 4 



C" 

= — — smd(x). 

Vp 



Thcse two wave functions must be consistent. This means that near the neighborhoods a, b of x 



Or equivalcntly, 



1 f b TT 

a(x) + 0(x) = -J pdx + - = (n+ 1)tt, n = 0, 1, 2, 3 • • • . 
(j) pdx = (n + n = 0, 1, 2, 3 ■ ■ ■ . 



For the half infinite potential well as shown in Fig. 4(b), one has 



pcía; = (n + n = 0, 1, 2, 3 • 



(72) 



(73) 



(74) 



(75) 



Analogously, the matching rule of the wave functions gives the quantization rule for the system with an infinitely deep 
square-well potential as illustrated in Fig. 4(c). Indeed, one has 



pdx = (n + í)h, n = 0, 1, 2, 3 • 



(76) 



The argument leading to the same result for a more general condition beyond the above examples can be found with 
the help of the Maslov index shown in the Rcf. [25]. 



(a) 



(b) 



V(x) 





(c) 



V(x) 



FIG. 4. WKB wave function matching boundary conditions for three cases of potentials. 
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